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Abstract 

The Fokker-Planck equation has been very useful for studying dy- 
namic behavior of stochastic differential equations driven by Gaussian 
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noises. However, there are both theoretical and empirical reasons 
to consider similar equations driven by strongly non-Gaussian noises. 
In particular, they yield strongly non-Gaussian anomalous diffusion 
which seems to be relevant in different domains of Physics. We there- 
fore derive in this paper a Fractional Fokker-Planck equation for the 
probability distribution of particles whose motion is governed by a 
nonlinear Langevin-type equation, which is driven by a Levy-stable 
noise rather than a Gaussian. We obtain in fact a general result for 
a Markovian forcing. We also discuss the existence and uniqueness of 
the solution of the Fractional Fokker-Planck equation. 

Running title: Fractional Fokker-Planck Equation and Levy noises 

Correspondence should be addressed to D. Schertzer (fax: +33 1 
44275259, e-mail: schertze@ccr.jussieu.fr). 
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1 Introduction and motivation 



The Fokker-Planck equation is one of the most celebrated equations 
in Physics, since it has been very useful for studying Q the dynamic 
behavior of stochastic differential equations driven by Gaussian noises. 
However, it turns out that many physical phenomena bring this frame- 
work into question. For instance, it has been argued that diffusion by 
geophysical turbulence P, [|l ^j, [j| |(| f?J corresponds, loosely speaking, 
to a series of sticking (pauses), when the particle is trapped by a co- 
herent structure, and (fast) flights, when the particle moves in the jet 
flow. 

Although there have been some attempts || to analyze and quan- 
tify this behavior with the help of the classical Fokker-Planck equa- 
tion, i.e. assuming finite moments of all orders, some laboratory ex- 
periments ]|, ||, ||] or numerical simulations of geostrophic turbulence 
H show that this phenomenology could be rather a consequence of 
the presence of heavy tails (i.e. power law fall-off) for the probability 
distribution and a strong anisotropy with a clearly preferred direction 
of diffusion. One can conclude |J] that if they are additive processes, 
the corresponding walks are Levy motions. 

Let us recall that indeed stable Levy motions L(t) generalize the 
Brownian motion B(t) in the sense that first they are also motions (e.g. 

[Tl[| ) whose increments AL(t, At) = L(t + At) — L(t) are station- 
ary (therefore AL has no statistical dependence on t) and independent 
for any non overlapping time lags At. Therefore, L(t) corresponds to 
the sum of independent, identically distributed Levy stable variables 
U, ll Hi. The second common property is that these incre- 



ments satisfy a "stability property" : up to a rescaling and recentring, 
the sum of different steps has the same probability distribution as 
one of the steps. Levy stable variables are precisely defined by this 
property. The stability property implies in both cases a property of 
attraction: under rather general conditions a renormalized sum of in- 
dependent identically distributed variables converge to a stable law. 
Furthermore, there are no other attractive laws. This explains why 
the stable property is so important. The attraction property corre- 
sponds to a broad generalization of the central limit theorem, with 
the important difference that whereas the classical theorem (Gaussian 
case) is satisfied with the condition that the variance is finite, the 
convergence towards a Levy law is obtained with the condition that 
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not only the variance of the summands Xi is infinite, but also that 
all their moments of order q equal to or larger than a critical order 
a (0 < a < 2) are infinite. This critical order a is called the Levy 
stability index and corresponds to the exponent of the power-law of 
probability distribution tails: 



any s » 1 : Pr(\AL\ > s) ^ s a <-> any q > a : E(\X\ q ) = oo 

(1) 

where Pr denotes the probability, E{.) is the mathematical expecta- 
tion and s is a given (large) nonnegative threshold. This statistical 
divergence of a Levy motion is due to jumps, whereas a Brownian 
motion is almost surely continuous. 

This index is the most important of the four parameters defining 
a Levy stable law. The second one is the 'skewness' (3 {—1 < (3 < 1) 
which defines the degree of asymmetry of the law, which is maximal for 
(3 = — 1 or (3 = +1, and the law is symmetric when (3 = 0. In spite of 
its name and some common properties, (3 nevertheless does not corre- 
spond to the classical skewness of a quasi-Gaussian law. The latter is 
indeed undefined for a stable Levy law due to the above mentioned sta- 
tistical divergences. The center 7 corresponds to the statistical mean 
when defined (i.e. a > 1) and/or to the median when symmetric (i.e. 
(3 = 0). The scale parameter D {D > 0) corresponds to a generaliza- 
tion of the variance of the Gaussian case. More precisely, as discussed 
below, it corresponds to the intensity scale of the cumulant of (pos- 
sibly non integer) order a. It yields an anomalous JlTj generalization 
of the classical Einstein relation: Var[X{t) — X(to)] = 2D(t — to), 
where Var(.) denotes the variance. Finally, let us emphasize that the 
Gaussian case corresponds to the limit case a = 2, which also implies 
(3 = 0, i.e. no asymmetry. 

Further comments are now in order on the relevance of Levy mo- 
tions in Physics. On the one hand, claims in favor of the relevance of 
Levy motions have been made on many physical phenomena ranging 
from subrecoil laser cooling [18, |H| to diffusion by flows in porous 



media [g0|, |2l]] , including finance fluctuations |22|) |2^] , see Ref . |24], |5[ 
for other examples. Many systems indeed display a phenomenology 
rather similar to that we reported above on geostrophic turbulence. 

On the other hand, important questions have been raised. In par- 
ticular, Ref. Iffil questioned the resulting infinite variance of the ad- 
verting field for porous media. Indeed, it turns out that recent esti- 



4 



mates [27] of the power-law of the probability distributions of the hy- 
draulic conductivity yields an exponent a ~ 3.5. The question of finite 
variance might apply to other examples, in particular for atmospheric 
turbulence where different studies [28] yield a critical exponent a ~ 7 
for the wind field. Therefore, in spite of their clear phenomenological 
interest, the relevance of pure Levy motions could be questioned. 

The main goal of this paper is to clarify and define a framework 
adequate for handling motions more general than pure Levy motions 
and which are nevertheless generated by the latter. We will do it by 



building upon a series of rather recent works [E9, 30, 32, 33, 34, 35 



[l?], |36j which show that the probability density of particles moving 
with a Levy motion satisfies a generalized Fokker-Planck equation 
involving fractional orders of differentiation. Indeed, it could be first 
argued in a "very formal and phenomenological" manner |^9| that a 
fractional power of the Laplacian yields an anomalous scaling for the 
corresponding diffusion. 

A Fractional Fokker-Planck equation was obtained in a less formal 
manner by |3^, |33j in the framework of the continuous time random 
walks (CTRW's) model of anomalous diffusion [31]. However, this 
method does not involve directly a stable Levy process, but a walk 
sharing some behavior common with the latter, without being equiv- 
alent to it. A different Fractional Fokker-Planck equation was intro- 
duced [35] with the help of a phenomenological and interesting trans- 
formation of the classical Fick law into a fractional Fick law. However, 
it is not clear that its solution corresponds to a (nonnegative) proba- 



bility distribution. A rather distinct approach was followed by |32|, 17 1 
since it starts with a linear Langevin-like equation with random forces 
which are exact stable Levy processes, which can be symmetric as well 
as asymmetric, and with no limitation on the possible values of the 
Levy index a. The fundamental mathematical tool which is used 
is the second characteristic (or cumulant generating) function of the 
motion defined by this Langevin-like equation. The particular case of 
symmetric processes correspond to what was previously inferred by 



[|29|, 30, 33, |35j]. However, it was shown that in the more general case 
of asymmetric processes, a new non-trivial 'advective-diffusive' term 
appears. This is confirmed with the help of a re-interpretation of the 
characteristic function of a Levy motion |36| . 

We already discussed that theoretically and empirically the non 
finiteness of the variance could be questioned. There are two more 
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general questions: the inhomogeneities of the medium, which are first 
emphasized for the introduction of the Levy motions, are finally re- 
duced to a (homogeneous) distribution of times when the particle is 
strongly kicked. As soon as this representation is granted, the medium 
(and its properties) does not intervene any longer. This is very re- 
strictive and for instance incompatible with the multifractality of the 
medium (or of the diffusion) when observed. The second reason is 
that the underlying processes are thought to be strongly nonlinear, 
whereas the transport is modeled with the help of a (stochastic) lin- 
ear equation. 

Both the successes and limitations of the previous results plead in 
favor of investigating a local and nonlinear modeling with the help of 
Levy motions. This is the reason that we investigate the properties a 
of nonlinear Langevin-like equation forced by a Levy stable motion. 



2 Statement of the problem 

Further to our above discussion, we consider the following nonlinear 
Langevin-like equation for a stochastic (real) quantity X(t) (e.g. lo- 
cation of a particle): 

dX(t) = m(X(t),t) dt + cr(X(t),t) dL (2) 

where the driving source is a Levy stable motion L(t) instead of Brow- 
nian motion B{t). The latter case corresponds to the basis of stochas- 
tic calculus (e.g. Ref. [[}?]]) and the corresponding differential equation 
is often called the Ito-Skorokhod equation. The extension to Levy sta- 
ble motion L{t) is rather natural and straightforward (e.g. ]3q| ) due 
to the common properties of L{t) and B(t) that we discussed in Sect. 
0, i.e. their infinitesimal increments are independent identically dis- 
tributed and furthermore stable. 

More precisely the Ito stochastic calculus corresponds to consider 
that the dL is, similarly to dB, a forward increment in time (it should 
be understood as dL(t, dt) = L(t + dt) — L{t)). This means that the 
value of X at time t is determined by events prior to the application 
of the stochastic force dL(t), which acts only from time t to t + dt. 

The Eq. || can also be understood under its integral form: 

X(t) = X(t ) + f m(X(t),t)dt + f o-(X(t),t) dL (3) 
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where the last term corresponds to a stochastic integration of a stochas- 
tic process. The integration of a stochastic process 3>(t) (in the case 
of Eq. ^: <I>(t) = a(X(t),t)) by the Levy motion L, is rather straight- 
forward in the case of step processes 



$(t) = /or t G (t n , t n+1 ), n = 0, 1, ... iV - 1 

- N-l 

I = £ ^^(Wl) " (4) 

n=0 

and this rather suggestive definition is naturally extended to functional 
spaces in which the step processes are dense. 

In order to establish local properties, for instance the time evolu- 
tion of the probability of the particles, we will use the differential form 
(Eq. |2[), whereas []32| , 17 1 rather used the integral form (EqQ) which 



becomes cumbersome in the nonlinear case and is in fact useful only 
to establish global properties (SectJ^). 

After having emphasized the similarities between L(t) and B(t), 
it is important to underline the non-trivial consequences due to the 
fact, contrary to the Gaussian case which has all its moments finite, 
Levy motions have a finite critical order of divergence of statistical 
moments (0 < a < 2). These include the fact that the mathematical 
techniques which could be used can be rather distinct. For instance, 
our derivation will rely on the use of the second characteristic function 
of the increments Sect.[| instead of probabilities of the increments as 
done usually for the derivation of the classical Fokker-Planck equa- 
tion. An obvious reason is that the former are relatively simple (see 
Sect.0), while the latter are not, with the only exception of the three 
following cases: a = 2, (3 = 0; a = 1, f3 = 0; a = 1/2, fi = 1. The fun- 
damental reason is that both the stability property and the divergence 
of moments are related to the presence of a cumulant of non integer 
order a. In relation to this problem, the convenient I? Hilbert struc- 
ture of Gaussian processes is reduced to a L a Banach structure for 
stable Levy processes. This is particularly important for the integral 
equation Eq.||, when defining functional spaces where step processes 
are dense. 

The linear case, which is the hitherto studied case, corresponds to: 
m(X(t),t) = m = Const., a(X(t),t) = a = Const. (5) 
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X(t) — X(to) is also a Levy motion which has the same Levy stability 
index a as its increments, but with a different center or trend and 
scale or amplitude. 

In the nonlinear case, m(X(t),t) and a(X(t),t) are (possibly non- 
linear) functions of X(t) and t, which satisfy certain regularity con- 
straints to be discussed later (Sect.|9|). They correspond to inhomo- 
geneities of the medium, which were ignored in the linear case. As 
a possibly important, but simple example, let us mention the Levy 
extension of the so-called geometric Brownian motion, which is rather 
ubiquitous and for instance is at the core of the Black-Scholes model 
for option pricing: m(X(t),t) = mX{t) and a(X(t), t) = o~X(t), where 
a is the 'votality' constant of the price X(t) of a given stock share. 

We will demonstrate that: 

Proposition 1 The transition probability density: 

Vt > to : p(x, t\x , t ) = Pr{X(t) = x\X(t ) = x ) (6) 

corresponding to the nonlinear stochastic differential equa- 
tion (Eq^), with a ^ 1 or (3 = 0, is solution of the follow- 
ing Fractional Fokker-Planck equation: . 

d d 

—p(x,t\x ,t ) = - — (rya[x,t) +m(x,t))p(x,t\x ,t ) 
-D[(-A) a / 2 (\a(x,t)\ a p(x,t\x ,t )) 
+p UJ (a)^(-A)^- i y 2 (\a(x,t)r 1 a(x,t)p(x,t\x ,t ))} (7) 

with the initial condition: 

p(x,to\x ,t ) = 5(x - Xq) (8) 

where 5{x — xq) is the degenerate Dirac measure in xq and 
uj(a) is defined by: 

Trot 

uj{a) = tan — (9) 
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and where the fractional powers of the Laplacian A will be discussed in 
Sect.||. Proposition [j] and Eq.0 are for scalar processes (i.e. A = 
and their extension to vector processes will be discussed and presented 
in Sect|| One may note that the fractional diffusive isotropic opera- 
tor — (— A) a / 2 applies via a fractional diffusivity \cr(x, t)\ a , whereas the 
advective-diffusive term corresponds to a conjugate action of a frac- 
tional diffusive term — (— A)( a_1 )/ 2 |cr(a;, and a convective term 
J^a(x,t) on the transition probability. 

This Fractional Fokker-Planck equation will be established with 
the help of the much more general proposition: 

Proposition 2 The inverse Fourier transform of the sec- 

ond characteristic function or cumulant generating func- 
tion of the increments of a Markov process X(t) generates 
by convolution the Fokker-Planck equation of evolution of 
its transition probability p(x, t\xo, to). 

We will demonstrate this proposition in a straightforward, yet rig- 
orous way. More precisely, we will establish the following: 

dp. . . f OK 

— {x,t\x ,to)= / dy—[x - y\y,t)p{y,t\x ,t ) (10) 

where K is the inverse Fourier transform of the cumulant generating 
function of the increments. The K arguments will become explicit in 
Sect. | 



Eq.10 not only holds for processes with stationary and indepen- 
dent increments, as in the linear case (Eq.||) but also for any Markov 
process, including those defined by the non-linear Langevin-like equa- 



tion (Eq.g with m ^ Const., a ^ Const.). As a consequence of Eq.10 
we will demonstrate the following: 

Proposition 3 When the increment's cumulant generat- 

ing function of a Markov process X(t) is defined by its ex- 
pansion in cumulants C n , its Fokker-Planck equation is : 

dp (— 1)" d n 

— (x,t\xQ,t ) = n , g-^[Cn(x,t)p(x,t\x 0l t )} (11) 

n£j 
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An obviously sufficient condition of convergence is obtained when 
the set J of the orders of differentiation n is finite. This is true in par- 
ticular for Gaussian forcing: J = {1, 2}. It corresponds to the classical 
Fokker-Planck equation. On the other hand, J = N would correspond 
to an analytic expansion of cumulants. In spite of its interest, we will 
not discuss the latter case in this paper, nor its relationship to the 
classical Kramers-Moyal expansion (e.g. |pQ| ). 

Below, we concentrate on the case of a finite, but non analytic 
expansion: J = {l,a} (non integer a, < a < 2), since it corre- 
sponds to the Levy extension (Sect.^ and yields Prop.[l| with the help 
of fractional derivatives, as discussed in Sect.||. 

3 The cumulant generating function 
of the increments 

The first and second (conditional) characteristic functions are respec- 
tively the moment generating function Zx(k,t — iol^o^o) and the 
cumulant generating function Kx(k, t — to\xo, to), associated with the 
transition probability p(x,t\xo,to) of a process X(t). These are de- 
fined by the Fourier transform of the latter, with k being the conjugate 
variable of x — xq : 

F\p(x,t\x ,t )] = Z x (k,t-t \x ,t ) (12) 
= exp(K x (k,t - t \x ,to)) (13) 
= E[exp(ik(X(t)-X(t ))\X(t ) = x ] (14) 

where E[-\-} denote the conditional mathematical expectation, F and 
F^ 1 respectively the Fourier-transform and its inverse: 

/oo 
dx exp(ikx)f(x) (15) 
-oo 

/°° dk 
— exp(-ikx)f(k) (16) 
-oo 27T 

The corresponding quantities for increments 5X(5t) = X(t + 5t) — 
X(t), corresponding to a given time lag 5t > 0, are defined in a similar 
way: 
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F[p(x + 5x,t + 5t\x,t)] 



6Z x (k,St\x,t) (17) 
exp(5K x (k,St\x,t)) (18) 
E[exp{ik{X(t + St) - X(t))\X(t) =(if) 



where k is the conjugate variable of 5x. The cumulants of the incre- 
ments C n are the coefficients of the Taylor expansion of SKx- 

(ik) n 

SK x (k,St\x,t) = 5tY ±J-C n (x,t) + o(St) (20) 

As already mentioned, the classical case corresponds to an analytic 
expansion of SK X , i.e. JCN, whereas we will be interested by a finite 
but non-analytic expansion J = {1, a} (non integer a, < a < 2). 

4 Processes with stationary and inde- 
pendent increments 

Let us first consider the simple sub-case of a process with station- 
ary and independent increments. It corresponds to C n (x,t) = C n = 
Const, in Eqs.|ll], |20] and as already discussed in Sect, [j], it includes 
the linear case (Eq.|5|) of the Langevin-like equation (Eq.Q). 

However, we believe that the following derivation is not only some- 
what pedagogical on the role of the characteristic functions for the 
nonlinear case, but also terser than derivations previously presented 
for the linear case. 

The stationarity of the increments implies that the transition prob- 
ability depends only on the time and space lags, i.e.: 



and similarly, the characteristic functions of the increments are no 
longer conditioned, for instance: 



p(x, t\x ,t ) = p(x -x ,t- t ) 




Z x (k,t- t \x ,t ) 
K x (k,t- t \x ,t ) 



Z x (k,t- t ) 
K x (k,t-t ) 




(23) 
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On the other hand, the independence of the increments implies that 
the transition probabilities satisfy a convolution (over any possible 
intermediate position y) for any given time lag 5t: 

VSt > : p(x-x ,t + 5t-t ) = J dy p(x-y,5t)p(y-x ,t-t ) (24) 

and the corresponding characteristic functions merely factor (resp. 
add). Therefore, we have: 



Z x (k,t + 5t-t )-Z x (k,t-t ) = Z x (k,t-t )(5Z x (k,5t)-l) (25) 
This in turn leads to: 

Z x (k,t+5t-t )-Z x (k,t-t ) = Z x (k,t-t )5K x (k,5t) + o{5t) (26) 
Its inverse Fourier transform yields: 



p(x,t+5t\x ,t )-p(x,t\x ,t ) = I dyF 1 [5K x (k, 5t)]p(y-x ,t-t )+o(5t) 

(27) 

This demonstrates (in the limit 5t — > 0) Prop]2| and Eq.lC, as well 
as Prop|| since Eq]27] corresponds, with the help of Eq.20, to: 

(— 1)" f r n \ 

p(x,t+5t\x ,t )-p(x,t\x ,t ) =St22 — [C n / dyS^2yp(y,t\x ,to)]+o(6t) 

naj n - J 

(28) 

where 6™ denotes the n th derivative of the Dirac function. Therefore, 
we obtain: 

Q (—l) n d n 

— p{x, t\x , t ) = — Cn ~Q^, p ( x ' * l^o, *o) (29) 

n£j 



which corresponds to the linear case of Eq.ll 
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5 More general Markov processes 



In the case of a Markov process which does not have stationary and 
independent increments, there is no longer a simple convolution equa- 
tion (Eq. ^4|) of the transition probabilities, nor a simple factorization 
of characteristic functions (Eq.^). However, the former satisfies a 
generalized convolution equation which corresponds to the Chapman- 
Kolmorogorov identity Eq] valid for any Markov process X(t): 



VSt > : p(x, t + 6t\x , t )) = J dy p(x, t + 6t\y, t)p(y , t\x , t ) (30) 

which indeed reduces to a mere convolution (Eq. ^4|) in the case of 
processes with stationary and independent increments. This identity 
can be written under the equivalent form: 



p(x, t + St\x , t ) = J dyj ^e~ ik y +SK ^ k ' St ^p(y,t\x ,t ) (31) 
Noting that we have: 

p(x,t\x ,t ) = / dy p(y,t\x ,to) / —e~ lky (32) 



2vr 



we obtain: 



p(x,t+5t\xo,t )-p{x,t\xo,to) =St J dyF 1 [5K x (k,5t\y,t)]p(y,t\x ,to)+o 

(33) 

In the limit 5t — ► 0, this corresponds to Prop. |2| and Eq. |l(J When 
J C N, it yields with the help of EqM 



neJ 



5p(x,t\x ,t ) = 5tJ2 J d y^x- y [ y n 7 Cn(y,t)p(y,t\x ,t )]+o(6t) 



(34) 

The limit 5t — > corresponds to Eq.ll and demonstrates Prop. |3] 
for a Markow process. 
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6 Extension to fractional orders 



In the two previous sections (Sects .Q- the fact that the indices 
n £ J should be integers intervene at best only in the correspondence 
between (integer order) differentiation -j^ (in Eq. |ll]) and powers of 
the conjugate variable k n (in Eq. |20| ). However, by the very definition 
of fractional differentiation (e.g. [f4l| ), this correspondence holds also 
for non integer orders. However, there is not a unique definition of 
fractional differentiation and therefore, as discussed in some details in 
|fT?j] ), we cannot expect to have a unique expression of the Fractional 
Fokker-Planck equation. 

Since in the following it will be sufficient to consider an expansion 
of the characteristic function involving fractional powers of only the 
wave number |fe|, it is interesting to consider Riesz's definition of a 
fractional differentiation. Indeed, the latter corresponds to consider 
fractional powers of the Laplacian: 

-(-Ar/ 2 f(x)=F- 1 [\k\ a f(k)] (35) 

which has furthermore the advantage of being valid for the vector 
cases. However, we will see in Sect. || that in general it does not apply 
in a straightforward manner for vector stable Levy motions. Indeed 
the latter introduces rather (one-dimensional) directional Laplacians, 
i.e. (one-dimensional) Laplacians along a given direction u (| u |= 1): 

- (-A^f(x) = F-^uTfik)] (36) 

where(.,.) denotes the scalar product. On the other hand, it will 
be useful to consider the fractional power of the contraction of the 
Laplacian tensor A: 



by a tensor a (a* denotes its transpose), with the following definition: 
- (-A : a.a*)f = (k,a.a*.k |f ] = a* .k \ a ] (38) 
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7 Levy case 



The second characteristic function of the increments 5L of the (scalar) 
Levy forcing is the following: 

k 

5K L (k, St) = &t[ihy - D\k\ a (l - iPjr-Mk, a)} + o(tft) (39) 

| A; | 

where u(k, a) is defined by: 

7TCt 7T 

a / 1 : uj(k, a) = u)(a) = tan—^-; a = 1 : u(k, a) = —log\k\ (40) 

Considering an Ito-like forward integration of Eq.||, the increments 
5L generates the following (first) characteristic function for the incre- 
ments 5X of the motion X(t): 



SZ x (k,St\x-Sx,t) = E{e ikm ^)8Z uL {k,8t\x,t) + +o{5t) (41) 

which yields the following elementary cumulant generating function 
5K X : " 



5Kx(k, 5t\x, t) = 5t[ikm(x, t) + ik , ya(x, t) 
-D\k\°\<T(x,t)ni - i(3 ^^ u;(k,a))]+o(St) (42) 

and which is of the same type as Eq.|20|, with J = {l,a}. Therefore, 
as discussed in Sectj6|, we have fractional differentiations in the corre- 
sponding Eqjll], which will precisely correspond to Eq.^, and therefore 
establishes Prop. |l]. 



8 Extension to vector processes 

With but one important exception, the extension of the previous re- 
sults to higher dimensions is rather straightforward. The starting 
point of this extension is the following nonlinear stochastic equation 
(X(t) G R d ): 
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dX(t) = m(X(t),t)dt + g(X(t),t).dL (43) 

where m is the natural d-dimensional vector extension of the deterministic- 
like trend, a is the dxd'-dimensional tensor extension of the modu- 
lation of the random driving force, and L is a d '-dimensional Levy 
stable motion. As discussed below, the expression of the characteris- 
tic function of the latter corresponds to the source of the difficulty in 
extending the scalar results to higher dimensions. On the contrary, 
it is straightforward to check that Props. |3| are valid in the vector 
case, with the following extensions (x G R d ) for Eq. [H]: 

^{x,t\x ,t ) = J dy^(x-y\y,t)p(y,t\x ,t ) (44) 
and for Eq. O (n G J C N d , | n \= £f =1 m): 



dp ^ _ x , (-l)lal f ^ f 

M { -> *'- ' t0 > ~ ^ (n 1 )\(n 2 )\..(n d )\dx?dx?..dx? t)P ^ %0 ' *° )J 

(45) 

On the other hand, Eq. 43 yields the following extension to Eq.f4l|: 

6Zx(k, 6t\x, t) = j <5Z g x(fe 5t\x, t) (46) 

and therefore we have: 



5Kx(k, 5t\x, t) = i k.m(x, t) + 5K L (g*.k, 5t\x, t) + o(5t) (47) 

Let us recall that a stable Levy vector in the classical sense [|l2| , f42| 
(see [44] for a discussion on a rather straightforward generalization, 



or |45|, 46, 47] for a more abstract generalization) corresponds to the 
limit of a sum of jumps, with a power-law distribution, along random 
directions u G dBi, B\ being the unit ball, distributed according to 
a (positive) measure dS(u). The latter, which generalizes the scale 
parameter D of the scalar case, is the source of the difficulty since in 
general the probability distribution of a stable Levy vector depends on 
this measure, and therefore is a non parametric distribution. However, 
as discussed below, there is at least a trivial exception: the case of 
isotropic stable Levy vectors. 
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Corresponding to our previous remarks, a (classical) stable Levy 
vector has the following (Fourier) cumulant generating function: 

K L (k) = 6t[i(Js.,i) - [ (ik,u) a dZ{u)} + o{6t) (48) 
which yields with the help of the Eq.|47|: 

d ~ r 

^K x (k) = -div(m + a.i)-F- 1 [ (iv*{x,t).k,u) a dZ(u)] (49) 

ot — JuedBi — 

The scalar case (Eq.|39[) corresponds to: 



7TCt 

0<p<l: (3 = 2p-l, dS(«) = Dcos( — )[p* (u _i) + (1 - p)5 (u+ i)] 

(50) 

For any dimension d, the second term on the right hand side of 



Eq.[49| corresponds to a fractional differentiation operator of order a. 
This operator can be slightly re-arranged. With the help of the odd 
dYi~(u) and even dT, + (u) parts of the measure dT,(u), 

2 dY, + (u) = dE(u) + dY l (-u); 2 dS~(u) = dS(u) - cffi(-u) (51) 
and the identity (0 being the Heaviside function): 

(ik) a = \k\ a [9(k)e iS f + 6{-k)e~ iS T] (52) 
one can write the extension of Eq.|7| under the following form: 

d 

— p(x,t\x ,t ) = -div[m(x,t) + g(x, t).-y)]p(x, t\x , t )] 

-[< (-A-g-g*) f >s+ - < (£■£*)■ (-4 -g-g*) ^ >^-]p(x,t\x o ,t )(53) 

where the symmetric fractional diffusive and respectively the antisym- 
metric advective-diffusive terms are defined , similarly to Eq. 35 , in the 
following manner: 

- < (-A :a.a*)% > s+ = / dZ+ivjF- 1 ^ {a* (x,t).k,u) \ a ] (54) 

- < (V.g*).(-4 : g-g*)^ >s- = 

/ dS'dO^IH^fetJ.&a) I (2* (£,*)•&, a) r -1 ] (55) 

Ju£dB 1 — 
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In general, each term corresponds to a rather complex integra- 
tion (which is indicated by the symbol < . >s)by the measure dE 



of directional fractional Laplacians (Eq.36). However, the symmetric 
term becomes simpler as soon as the even part dS + of the measure 
<i£ is isotropic. Indeed, the integration over directions yields only a 
prefactor D: 



< -(A : > s += D (-A : a.a*)' 



D = / dZ + (u) | (u^u) \ a (56) 

and for a = 2 this corresponds to the classical term (A : a. a*) of the 
standard d-dimensional Fokker-Planck equation. If dY> itself is rota- 
tion invariant, then the asymmetric operator vanishes, since dS~ = 0. 
If furthermore, a is scalar, i.e. a = al, then one obtains the following 
Fractional Fokker-Planck equation: 



d 

— p{x,t\x ,t ) = -div\a.2(x,t) + m(x,t)]p(x,t\x ,t ) (57) 

- D [(-A) a / 2 ]|a(s,t)| a p(£,t|so,to) (58) 

Therefore, as one might expect it, the rotation symmetries yields 
a rather trivial extension of the standard Gaussian case: a fractional 
power a of the d-dimensional Laplacian, as in the pure scalar case 
(Eq.0). Obviously, the integration performed in Eq.53 is also greatly 



simplified as soon as dH(u) is discrete, i.e. its support corresponds to 
a discrete set of directions . 

On the other hand, let us note that the framework of generalized 



stable Levy vectors |44|, [46|, |47[], allows one to introduce a much 



stronger anisotropy than the the measure cE does it for classical stable 
Levy vectors. This therefore diminishes the importance of the asym- 
metry of the latter. Indeed, the components of a generalized stable 
Levy vector do not necessarily have the same Levy stability index, 
the latter being generalized into a second rank tensor. Similarly, the 
differential operators involved in the corresponding Fractional Fokker- 
Planck equation no longer have a unique order of differentiation. This 
is rather easy to check in case of a discrete measure dS (u) and we will 
explore the general case elsewhere. 
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9 Existence and uniqueness of the so- 
lution 



The previous sections established a generalization of the Fokker-Planck 
equation for the evolution of the probability distribution of nonlinear 
stochastic differential equations driven by Levy stable noises. This 
was the main goal of this paper. Naturally, one would also like to 
have if possible a theorem of existence and uniqueness of the solution 
of this equation. Due to its origin, such a theorem will also imply 
that the solution will remain positive and normalized, as required for 
a transition probability. In this section we argue that the general 
results obtained jl8| in the classical Gaussian case (a = 2) are also 
relevant for the Levy extension, whereas up until now existence and 
uniqueness conditions of partial fractional differential equations have 
been scarcely explored (see however Ref . |^9| , Henry) and therefore we 
cannot rely on general results. 

The classical Fokker-Planck equation belongs to the well explored 
domain of parabolic equations. Existence and uniqueness of the solu- 
tion fundamentally result fl5l| from the fact that the linear operator 
A = —A is a (self- adjoint) positive generator of a semigroup of con- 
traction operators T(t) = e~ tA , t > 0. In the case of constant coef- 
ficients (linear Langevin equation), the solution is directly obtained 
with the help of T(t) and this ensures its existence, uniqueness and 
positiveness. Note that in our case, the semigroup action corresponds 
to the equation of convolution (Eq,p4|). 

Similar results hold for a Lipschitz variation of the coefficients, i.e.: 

\m(x, t) — m(y, t)\ + \a(x, t) — a(y, t)\ < D\x — y\ (59) 

as well as a condition of slow growth in time of the coefficients m(X(t),t) 
and a(X(t), t),e.g.: 

\m(x,t)\ + \a(x,t)\ < C|l + x\ (60) 

where D and C are given positive constants. 

These conditions have been extensively used for the classical Fokker- 
Planck equation with non constant coefficients (e.g. ^0|). Consider- 
ing now the fractional generalization, it is important to note that the 
fractional power of the Laplacian — (— A) a / 2 remains positive, since its 
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definition Eq.35 corresponds to replacing the eigenvalues k 2 by eigen- 
values having as real part \k\ a . Therefore, we remain inside of the 
previous framework of contraction semigroup and the previous results 
should hold. 

This could be also seen on the integral form of the differential 
equation. Indeed, in the classical case, the Lipschitz condition is clas- 
sical for the Brownian forcing [52, 38], as well as for the more general 
case of martingale and semi-martingale forcing [53, 54, 55]. The latest 
case is relevant for the stable Levy forcing. The Lipschitz condition 
can be rather understood as a condition of convergence of the Picard 
iteration method towards a fixed point: 



X n+1 {t) =X(t )+ j m(X n (t),t)dt + j a{X n (t),t) dL; X° '(t) = X (t ) 

(61) 

On the other hand, the condition of slow growth (^0|) in time pre- 
vents a finite explosion time, i.e. X(t) remains finite for any given 
finite time t: this condition is rather general, since it is already re- 
quired by the deterministic part of the Langevin-like equation. 



10 Conclusion 

We have derived a Fractional Fokker-Planck equation, i.e. a kinetic 
equation which involves fractional derivatives, for the evolution of the 
probability distribution of nonlinear stochastic differential equations 
driven by non-Gaussian Levy stable noises. We first established this 
equation in the scalar case, where it has a rather compact expression 
with the help of fractional powers of the Laplacian, and then dis- 
cussed and presented its extension to the vector case. This Fractional 
Fokker-Planck equation generalizes broadly previous results obtained 
for a linear Langevin-like equation with a Levy forcing, as well as the 
standard Fokker-Planck equation for a nonlinear Langevin equation 
with a Gaussian forcing. 
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